1 



STAR FORMATION: 3D COLLAPSE OF TURBULENT CLOUD CORES 
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Abstract 

We have performed fully 3-D simulations of the col- 
lapse of molecular cloud cores which obey the logatropic 
equation of state. By following the collapse of these cores 
from states of near hydrodynamic equilibrium, we are 
able, to produce accretion histories which closely resem 



degree of turbulent support via a phenomenological equa- 
tion of state. This logatropic equation of state reads 
P/P c = 1 + A hx(p/p c ) where the subscript "c" de- 
notes "central" values. The constant A is an adjustable 
parameter whose value for a range of molecular cloud 



blc tlhose of observer! cotes The accretion proceeds in fm.ir fc Pudritz (1996)| to be A 



core masses was determined empirically by McLaughlin 

0.2 ±0.02. 



distinct stages: an initial period of very slow accretion; a 
period of vigorous accretion following the development of 
a central density singularity, with M « t 3 , as predicted by 
self-similar models; a period of relatively stable, vigorous 
accretion; and finally a gradual decrease in the accretion 
rate once about 50% of the available mass of the molec- 
ular cloud core has been accreted. These results may ex- 
plain the accretion histories of cores as they pass through 
the pre-protostellar, Class 0, and Class I stages. 
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1. Introduction 

Mounting evidence indicates that massive molecular cloud 
cores possess properties which are not adequately ex- 
plained by models employing an isothermal equation of 
state. For example, surveys of prestellar cores in a va- 
riety of star-forming complexes show that intermediate 
and high mass cores have nonthermal line widths which 
are significantly greater than their thermal line widths 



(Caselli & Myers 1995). Models of the collapse of isother- 



mal molecular cloud cores exhibit a number of features 
which are not borne out by observations, such as a t ime 
invar iant rate of accretion onto the central object (Shu 



1993| ) 



1977 ) and highly supersonic infall ( Foster fe Chevalier 



To address some of the shortcomings of isothermal 



The logatropic equation of state admits two equilib- 
rium solutions. The features of each solution are devel- 
oped in McLaughlin & Pudritz (1996). The "singular" 



models, McLaughlin & Pudritz 1996 advanced a model 



equilibrium solution is analogous to the singular isother- 
mal sphere (SIS) of 3hu (1977), except that its density 
profile is p oc r . The "nonsingular" solution is analo- 
gous to the Bonnor-Ebert sphere, in that it is a pressure- 
truncated equilibrium solution with a finite central den- 
sity and density profile in its outer regions which matches 
that of the singular solution. 

The nonsingular logatrope has the potential to be 
as useful a model for massive star-forming cores as the 
Bonnor-Ebert sphere is for low mass star-forming cores 
(see, for example, Alves, Lada, fc Lada (2001) , who showed 
that a Bonnor-Ebert sphere is an excellent fit to the struc- 
ture of dark cloud B68). Hence, we undertook a numerical 
study of the collapse of a nonsingular logatropic sphere 
from an initially near-equilibrium state until most of its 
mass had been accreted onto the central condensation. 



2. Numerical Methods and Tests 

Our fully 3D simulations were conducted using the ZEUS- 
MP parallelized hydrodynamics code. We used a 256 3 
Cartesian grid with a "sink cell" at its centre. Mass falling 
into the sink cell was presumed to have accreted onto 
the central protostar and was represented thereafter by a 
point mass at the grid centre. The initial configuration for 
each of our simulations was a spherical molecular cloud 
core in near hydrodynamic equilibrium, surrounded by an 
external medium of low and constant pressure. We con- 
ducted extensive tests of our numerical method using the 



for molecular cloud core collapse which accounts for some ;ollapse of the singular logatrope (McLaughlin & Pudritz 



Proceedings of the conference "Chemistry as a Diagnostic of Star Formation, " University of Waterloo, Canada, 21-23 
August 2002 (C.L. Curry & M. Fich eds.). 



2 



Reid, Pudritz, and Wadsley: 3D Collapse of Turbulent Cloud Cores 



1997) as our test case. Excellent agreement betwe en the 
analytic and numerical solutions was found, (see Reid. 
Pud ritz, fc Wadsley (2002) for a full description of our 
numerical method). 



3. Results 

The most interesting of our results from the collapse of 
the nonsingular logatrope are presented in Figures |^-|[ 
Figures [l] and || show the evolution of the density and 
infall speed profiles of a collapsing nonsingular logatrope 
prior to the formation of a singular density profile. The 
model has been scaled such that the total mass of the core 
is 1 Mq . The total core mass is determined by the choice 
of central temperature and external pressure — our choices 
of T c — 10 K and P s = 1.3x 10 5 result in a core of total 
mass 1 M Q , but the results can be rescaled and applied to 
models of higher mass (for example, by assuming a higher 
confining pressure). Indeed, observational evidence has 
emerged which indicates that the confining pressures on 
molecular cloud cores are one or two orders of magnitude 
higher than we assumed ( Johnstone ct al. 2000] ). 

Perhaps one of the most notable features of Figure 
[j] is the behaviour of the density profile at late times. 
At the moment when the density becomes singular (la- 
beled t = 0), the density in the inner parts of the core 
scales as r~ 3 / 2 . This is the density profile expected for 
a collapsing singular logatropic sphere, as well as for a 
collapsing SIS. This implies that, whereas observations 
of r~ 3 / 2 density profiles in molecular cloud cores have 
traditionally been interpreted as evidence for isothermal 
collapse, this should no longer be the case. Even mod- 
els which incorporate turbulent support will show such a 
density profile. 

Because of their extra nonthermal support, nonsin- 
gular logatropic spheres collapse more "gently" than do 
their isothermal counterparts. As can be seen by com- 
paring our Figure with Figure la of Foster & Chevalier 



(1993), the logatropic collapse is more gentle than that 
of the Bonnor-Ebert sphere. Foster & Chevalier (1993) 
find that, at the time of singularity formation, 44% of 
the mass of a Bonnor-Ebert sphere is in supersonic infall. 
We find that less than 5.5% of a nonsingular logatrope is 
in supersonic infall at the same point in its evolution. 

Figure || shows the rate of accretion onto the central 
object in collapsing nonsingular logatropes of masses 1, 
2.5, and 5 M . The reader should bear in mind that these 
are fiducial masses only — the models can be rescaled to 
produce much more massive cores and hence higher ac- 
cretion rates by increasing the confining pressure. Each 
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Figure 1. Density evolution of the 1 Mq nonsingular core. 
The singular density profile is reached at t=0. The thick 
straight line indicates p oc r" 1 and the thinner straight 
line indicates p oc r~ 3 / 2 . The solid curve (t = —1.094^ 
represents the initial density profile of the nonsingular 
core. Adapted from Reid, Pudritz, & Wadsley (2002). 

curve shows the evolution of the accretion rate from the 
near-equilibrium starting point and until 95%, 91%, and 
29% of the total available mass has been accreted (for the 
1, 2.5, and 5 M Q cases respectively). 

4. Discussion 

The accretion profiles shown in Figure || bear several 
important similarities to what is known observationally 
about accretion in young stellar objects. A simple pic- 
ture of the evolution of an accreting molecular cloud core 
suggests an accretion history of the following form: accre- 
tion proceeds slowly through the pre-protostellar (PPC) 
phase, increases markedly at the formation of a central 
source (to power the substantial outflows seen in Class O 
sources), and finally declines and dissipates as the core 
passes into the Class I stage. Andre, Ward- Thompson, & 



Barsony (2000) mark the transition between the Class O 



and Class I phase at the point where the mass accreted 
onto the central object is roughly equal to that remaining 
in the envelope. Our results show a similar transition at 
the same point. As indicated for the 1 M Q curve in Figure 
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Figure 2. Evolution of the radial infall speed, \v r \, of a 
1 Mq nonsingular logatrope preceding the development of 
the singular density profile. Adapted from Reid, Pudritz, 
& Wadsley (2002). 



there is an initial period of slow accretion, which we 
liken to the PPC stage. As soon as the density profile of 
the core becomes singular (t = 0), which can be thought 
of as the moment of protostar formation, the accretion 
rate increases markedly, launching a phase of vigorous ac- 
cretion sufficient to power an outflow. We liken this to the 
Class phase. Once half of the available mass has been 
accreted, the accretion rate declines steadily to near-zero. 

We suggest that our simulations of the collapse of lo- 
gatropic molecular cloud cores bear significant similarities 
to the observations. They offer an alternative method for 
modelling the structure of intermediate and higher mass 
molecular cloud cores. Cores which exhibit density power- 
law indices of between 1 and 1.5 (see, for example, van 
der Tak et ah (2000)) should be pursued as possible can- 
didate nonsingular logatropes. Our results bolster the ar- 
gument that models of collapsing protostellar cores must 
account for some degree of nonthermal support in order 
to reproduce the observations. 
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FigureS. Time evolution of M acc (t) and M acc (t) for the 
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